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A SYSTEM AND METHOD FOR COMPUTING LOW COMPLEXITY 
ALGEBRAIC NETWORK CODES FOR A MULTICAST NETWORK 

10 BACKGROUND 
Technical Field: 

The invention is related to automatically optimizing network coding, and in 
15 particular, to a system and method for performing network coding at nodes in the 
network along transmission paths between senders and receivers for increasing 
an amount of information that can be reliably broadcast from a sender to a 
collection of receivers through the network. 

20 Related Art: 

In general, a multicast network can be described as a network with 
directed edges. There are a number of existing schemes for routing network 
flows in an attempt to optimize the capacity of such networks. 

25 

For example, one conventional scheme has demonstrated that if coding is 
allowed at internal nodes in a network, then the multicast capacity is in general 
higher than if no such coding is allowed. Further, this scheme has also 
demonstrated the existence of multicast codes that would achieve a natural 
30 upper bound on a multicast capacity by applying a max-flow min-cut theorem to 
the network between a sender and a number of receivers. Unfortunately, the 
results offered by this scheme depend on random coding arguments without 
providing any construction techniques for practical multicast codes. 
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Another fairly recently offered scheme for coding acyclic networks has 
demonstrated a connection between algebraic geometry and network coding. 
This network coding scheme examined the performance of codes where nodes 
are allowed to group together incoming bits into blocks of a predetermined 

5 length, m. The resulting symbols are then treated as elements in a finite field 
having a size of 2 m . This scheme then performs a linear combination on the 
symbols in the finite field to produce outgoing symbols which are elements in a 
finite field. Decoding at receiver nodes is also a linear operation over the finite 
field on the incoming symbols. This scheme also provides techniques for 

10 examining multicast scenarios, such as, for example, certain edge failure 

patterns may occur, networks with delay, and other special encoding scenarios 
(such as when all sources wish to transmit all their information to all sinks in the 
network). 

15 However, this scheme for coding acyclic networks has several drawbacks. 

For example, where N represents a number of receivers in the network and C 
represents a cutset capacity of the network, the solution to the multicast problem 
requires a field size q for the network codes to be larger then NC. This number 
quickly becomes impractically large for arithmetic implementation as the size of N 

20 and C increase. Further, the resulting codes would involve "flooding" the 
network, thereby likely using more network edges than would otherwise be 
necessary for the same or greater network capacity. 

Therefore, what is needed is a system and method for coding networks 
25 that overcomes the disadvantages of the aforementioned schemes. For 

example, such a system and method should provide construction techniques for 
practical multicast codes. Further, these multicast code construction techniques 
should limit the complexity of network coding, even on large networks, such that 
any arithmetic computations are feasible. Finally, the resulting network codes 
30 should avoid flooding of the network in order to optimize capacity. 
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SUMMARY 



A "multicast network code constructor," or simply "network code 
constructor," as described herein, provides a system and method for facilitating 

5 network based coding in a multicast environment by constructing efficient 

multicast network codes for optimizing network flows, thereby increasing reliable 
network throughput. In general, the network codes are used to process incoming 
data at each node on a byte-by-byte level to produce outgoing packets to each 
neighboring node in the multicast network. The network code constructor 

10 provides for multicast network coding, or simply "network coding," in which 

arithmetic operations can occur in any finite field with more than N-l elements, 
where N represents the number of receivers in the network. Further, the 
complexity of arithmetic employed by the network code constructor is 
independent of the network capacity, and dependent only on the number of 

15 receivers in the network. In addition, in one embodiment, multicast network 
codes are restricted to the portion of the network obtained by a union of unicast 
flows from a sender node to each receiver node to produce codes which do not 
flood the network excessively, and thus have an associated lower design 
complexity. 

20 

The multicast networks described herein are understood to be defined as 
networks with directed edges that are assumed to be error-free, with a single 
sender of information s wishing to transmit the same information to N receivers 
K = {r v r 2i ...,r N } . However, it should be noted that while the descriptions of the 
25 network code constructor provided herein are explained in the context of an 
analysis of multicast networks over directed acyclic graphs with zero-delay 
edges, the ideas described herein are easily extensible to robust network codes, 
more general encoding scenarios, and networks with delays. 

30 The network code constructor provides an efficient coding process in 

which arithmetic operations can occur in any finite field with more than N-l 
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elements. Thus the complexity of arithmetic is independent of the capacity of the 
network, and dependent only on the number of receivers. The coding design 
process described herein provides a solution to a multicast network coding 
problem, wherein the complexity of the solution is dominated by the time 
5 complexity of running a maximum flow algorithm N times, plus 

0(N(N + C) 3 1 E |) operations over any finite field with more than N - 1 elements, 
where E represents the number of edges. 

In general, the network code constructor can be simply described as a 
10 system and method for computing a linear network code. Computing the network 
code is accomplished by first computing flows between at least one sender and 
two or more receivers given capacities on each network edge. Next, given the 
computed flows, network node coefficients are then computed directly from the 
computed flows to provide an efficient and computationally inexpensive network 
15 coding, relative to conventional network coding schemes. 

In particular, the network code constructor facilitates network based 
coding in a multicast environment by determining efficient codes for optimizing 
network flows. This process begins by using known parameters of the network, 

20 including, for example, information describing a sender, S, receivers, n and each 
internal node in the network, along with the flow capacity (edges) of each of 
these network elements. These parameters are then used to reduce the network 
to a form that is more conducive to analysis. In particular, the network is reduced 
to a network whose edges have unit capacities by replacing each edge having a 

25 capacity c with c edges having unit capacity. 

Following this reduction of the network, a conventional max-flow algorithm 
is used to determine the multicast capacity C from the sender to the receivers. 
Then, given this multicast capacity C, along with the known network parameters, 
30 the network code constructor constructs multicast network codes by computing a 
flow of rate R < C from the sender to each receiver, and from these flows directly 
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computes encoding vectors and decoding matrices for transmitting symbols from 
the sender to the receivers across the network, and then decoding those symbols 
at each receiver. 

5 In addition to the just described benefits, other advantages of the system 

and method for automatically determining efficient codes for optimizing network 
flows will become apparent from the detailed description which follows 
hereinafter when taken in conjunction with the accompanying drawing figures. 

10 

DESCRIPTION OF THE DRAWINGS 

The specific features, aspects, and advantages of the present invention 
will become better understood with regard to the following description, appended 
15 claims, and accompanying drawings where: 

FIG. 1 is a general system diagram depicting a general-purpose 
computing device constituting an exemplary system for automatically facilitating 
network based coding in a multicast environment by determining efficient codes 
20 for optimizing network flows. 

FIG. 2 illustrates an exemplary architectural diagram showing exemplary 
program modules for automatically facilitating network based coding in a 
multicast environment by determining efficient codes for optimizing network 
25 flows. 

FIG. 3 illustrates an exemplary system flow diagram for automatically 
reducing a network with integer edge capacities to a network with unit edge 
capacities. 

30 
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FIG. 4 illustrates an exemplary system flow diagram for automatically 
determining a multicast capacity of a network. 

FIG. 5 illustrates an exemplary system flow diagram of an initialization 
5 stage for automatically designing a code for a network. 

FIG. 6 illustrates an exemplary system flow diagram of an encoder design 
stage for automatically designing a code for a network. 

10 FIG. 7 illustrates an exemplary system flow diagram of a decoder design 

stage for automatically designing a code for a network. 

FIG. 8 illustrates an exemplary system flow diagram for automatically 
computing encoding vectors for every edge of a multicast network. 

15 

FIG. 9 illustrates an exemplary system flow diagram for automatically 
choosing a linear combination of vectors for use in a greedy algorithm utilized by 
the network code constructor. 

20 FIG. 10 illustrates an exemplary system flow diagram for automatically 

multicasting symbols from a source to every receiver using network coding 
generated by the network code constructor, as illustrated in FIG. 3 through FIG. 
9. 

25 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 

In the following description of the preferred embodiments of the present 
invention, reference is made to the accompanying drawings, which form a part 
30 hereof, and in which is shown by way of illustration specific embodiments in 
which the invention may be practiced. It is understood that other embodiments 
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may be utilized and structural changes may be made without departing from the 
scope of the present invention. 

1.0 Exemplary Operating Environment: 

5 

Figure 1 illustrates an example of a suitable computing system 
environment 100 on which the invention may be implemented. The computing 
system environment 100 is only one example of a suitable computing 
environment and is not intended to suggest any limitation as to the scope of use 
10 or functionality of the invention. Neither should the computing environment 100 
be interpreted as having any dependency or requirement relating to any one or 
combination of components illustrated in the exemplary operating environment 
100. 

15 The invention is operational with numerous other general purpose or 

special purpose computing system environments or configurations. Examples of 
well known computing systems, environments, and/or configurations that may be 
suitable for use with the invention include, but are not limited to, personal 
computers, server computers, hand-held, laptop or mobile computer or 

20 communications devices such as cell phones and PDA's, multiprocessor 
systems, microprocessor-based systems, set top boxes, programmable 
consumer electronics, network PCs, minicomputers, mainframe computers, 
distributed computing environments that include any of the above systems or 
devices, and the like. 

25 

The invention may be described in the general context of computer- 
executable instructions, such as program modules, being executed by a 
computer. Generally, program modules include routines, programs, objects, 
components, data structures, etc., that perform particular tasks or implement 
30 particular abstract data types. The invention may also be practiced in distributed 
computing environments where tasks are performed by remote processing 



-7- 



devices that are linked through a communications network. In a distributed 
computing environment, program modules may be located in both local and 
remote computer storage media including memory storage devices. With 
reference to Figure 1 , an exemplary system for implementing the invention 
5 includes a general-purpose computing device in the form of a computer 110. 

Components of computer 110 may include, but are not limited to, a 
processing unit 120, a system memory 130, and a system bus 121 that couples 
various system components including the system memory to the processing unit 

10 120. The system bus 121 may be any of several types of bus structures 

including a memory bus or memory controller, a peripheral bus, and a local bus 
using any of a variety of bus architectures. By way of example, and not 
limitation, such architectures include Industry Standard Architecture (ISA) bus, 
Micro Channel Architecture (MCA) bus, Enhanced ISA (EISA) bus, Video 

15 Electronics Standards Association (VESA) local bus, and Peripheral Component 
Interconnect (PCI) bus also known as Mezzanine bus. 

Computer 110 typically includes a variety of computer readable media. 
Computer readable media can be any available media that can be accessed by 

20 computer 110 and includes both volatile and nonvolatile media, removable and 
non-removable media. By way of example, and not limitation, computer readable 
media may comprise computer storage media and communication media. 
Computer storage media includes volatile and nonvolatile removable and non- 
removable media implemented in any method or technology for storage of 

25 information such as computer readable instructions, data structures, program 
modules, or other data. 

Computer storage media includes, but is not limited to, RAM, ROM, 
EEPROM, flash memory, or other memory technology; CD-ROM, digital versatile 
30 disks (DVD), or other optical disk storage; magnetic cassettes, magnetic tape, 
magnetic disk storage, or other magnetic storage devices; or any other medium 



-8- 



which can be used to store the desired information and which can be accessed 
by computer 110. Communication media typically embodies computer readable 
instructions, data structures, program modules or other data in a modulated data 
signal such as a carrier wave or other transport mechanism and includes any 

5 information delivery media. The term "modulated data signal" means a signal 
that has one or more of its characteristics set or changed in such a manner as to 
encode information in the signal. By way of example, and not limitation, 
communication media includes wired media such as a wired network or direct- 
wired connection, and wireless media such as acoustic, RF, infrared, and other 

10 wireless media. Combinations of any of the above should also be included within 
the scope of computer readable media. 

The system memory 130 includes computer storage media in the form of 
volatile and/or nonvolatile memory such as read only memory (ROM) 131 and 

15 random access memory (RAM) 132. A basic input/output system 133 (BIOS), 
containing the basic routines that help to transfer information between elements 
within computer 110, such as during start-up, is typically stored in ROM 131. 
RAM 132 typically contains data and/or program modules that are immediately 
accessible to and/or presently being operated on by processing unit 120. By way 

20 of example, and not limitation, Figure 1 illustrates operating system 134, 

application programs 135, other program modules 136, and program data 137. 

The computer 110 may also include other removable/non-removable, 
volatile/nonvolatile computer storage media. By way of example only, Figure 1 

25 illustrates a hard disk drive 141 that reads from or writes to non-removable, 
nonvolatile magnetic media, a magnetic disk drive 151 that reads from or writes 
to a removable, nonvolatile magnetic disk 152, and an optical disk drive 155 that 
reads from or writes to a removable, nonvolatile optical disk 156 such as a CD 
ROM or other optical media. Other removable/non-removable, 

30 volatile/nonvolatile computer storage media that can be used in the exemplary 
operating environment include, but are not limited to, magnetic tape cassettes, 
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flash memory cards, digital versatile disks, digital video tape, solid state RAM, 
solid state ROM, and the like. The hard disk drive 141 is typically connected to 
the system bus 121 through a non-removable memory interface such as interface 
140, and magnetic disk drive 151 and optical disk drive 155 are typically 
5 connected to the system bus 121 by a removable memory interface, such as 
interface 150. 

The drives and their associated computer storage media discussed above 
and illustrated in Figure 1, provide storage of computer readable instructions, 

10 data structures, program modules and other data for the computer 110. In Figure 
1 , for example, hard disk drive 141 is illustrated as storing operating system 144, 
application programs 145, other program modules 146, and program data 147. 
Note that these components can either be the same as or different from 
operating system 134, application programs 135, other program modules 136, 

15 and program data 137. Operating system 144, application programs 145, other 
program modules 146, and program data 147 are given different numbers here to 
illustrate that, at a minimum, they are different copies. A user may enter 
commands and information into the computer 110 through input devices such as 
a keyboard 162 and pointing device 161, commonly referred to as a mouse, 

20 trackball, or touch pad. 

Other input devices (not shown) may include a microphone, joystick, game 
pad, satellite dish, scanner, radio receiver, and a television or broadcast video 
receiver, or the like. These and other input devices are often connected to the 

25 processing unit 120 through a user input interface 160 that is coupled to the 
system bus 121 , but may be connected by other interface and bus structures, 
such as, for example, a parallel port, game port, or a universal serial bus (USB). 
A monitor 191 or other type of display device is also connected to the system bus 
121 via an interface, such as a video interface 190. In addition to the monitor, 

30 computers may also include other peripheral output devices such as speakers 
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197 and printer 196, which may be connected through an output peripheral 
interface 1 95. 

The computer 1 10 may operate in a networked environment using logical 
5 connections to one or more remote computers, such as a remote computer 180. 
The remote computer 180 may be a personal computer, a server, a router, a 
network PC, a peer device, or other common network node, and typically 
includes many or all of the elements described above relative to the computer 
110, although only a memory storage device 181 has been illustrated in Figure 1. 
10 The logical connections depicted in Figure 1 include a local area network (LAN) 
171 and a wide area network (WAN) 173, but may also include other networks. 
Such networking environments are commonplace in offices, enterprise-wide 
computer networks, intranets, and the Internet. 

15 When used in a LAN networking environment, the computer 1 10 is 

connected to the LAN 171 through a network interface or adapter 170. When 
used in a WAN networking environment, the computer 110 typically includes a 
modem 172 or other means for establishing communications over the WAN 173, 
such as the Internet. The modem 172, which may be internal or external, may be 

20 connected to the system bus 121 via the user input interface 160, or other 
appropriate mechanism. In a networked environment, program modules 
depicted relative to the computer 1 10, or portions thereof, may be stored in the 
remote memory storage device. By way of example, and not limitation, Figure 1 
illustrates remote application programs 185 as residing on memory device 181. 

25 It will be appreciated that the network connections shown are exemplary and 
other means of establishing a communications link between the computers may 
be used. 

The exemplary operating environment having now been discussed, the 
30 remaining part of this description will be devoted to a discussion of the program 
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modules and processes embodying a system and method for automatically 
determining efficient codes for optimizing network flows. 

2.0 Introduction: 

5 

A "multicast network code constructor," or simply a "network code 
constructor," as described herein, provides a system and method for facilitating 
network based coding in a multicast environment by computing efficient multicast 
network codes for optimizing network flows, thereby increasing reliable network 

10 throughput. In general, the multicast network codes are used to process 

incoming data at each node on a byte-by-byte level to produce outgoing packets 
to each neighboring node in the multicast network. The network code 
constructor provides for multicast network coding, or simply "network coding," in 
which arithmetic operations can occur in any finite field with more than N-\ 

15 elements, where N represents the number of receivers in the network. Further, 
the complexity of arithmetic employed by the multicast network code constructor 
is independent of the network capacity, and dependent only on the number of 
receivers in the network. In addition, in one embodiment, multicast network 
codes are restricted to the portion of the network obtained by a union of unicast 

20 flows from a sender node to each receiver node to produce codes which do not 
flood the network excessively, and thus have an associated lower design 
complexity. 

The multicast networks described herein are understood to be defined as 
25 networks with directed edges that are assumed to be error-free, with a single 
sender of information s wishing to transmit the same information to N receivers 
TZ = {r v r 2v ..,r N } . However, it should be noted that while the descriptions of the 

network code constructor provided herein are explained in the context of an 
analysis of multicast networks over directed acyclic graphs with zero-delay 
30 edges, the ideas described herein are easily extensible to robust network codes, 
more general encoding scenarios, and networks with delays. 
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The network code constructor provides an efficient coding process in 
which arithmetic operations can occur in any finite field with more than N-l 
elements. Thus the complexity of arithmetic is independent of the capacity of the 
network, and dependent only on the number of receivers. The coding design 
5 process described herein provides a solution to a multicast network coding 
problem, wherein the complexity of the solution is dominated by the time 
complexity of running a maximum flow algorithm N times, plus 
0(N(N + C) 3 1 E |) operations over any finite field with more than N - 1 elements, 
where E represents the number of edges. 

10 

In general, the network code constructor can be simply described as a 
system and method for computing a linear network code. Computing the 
multicast network code is accomplished by first computing flows between at least 
one sender and two or more receivers given capacities on each network edge. 
15 Next, given the computed flows, network node coefficients are then computed 
directly from the computed flows to provide an efficient and computationally 
inexpensive network coding, relative to conventional network coding schemes. 

2.1 System Overview: 

20 

In general, the network code constructor operates by computing network 
codes for a known network, having senders and receivers of known flow 
capacity. The following section briefly describes the general operation of the 
multicast network code constructor in terms of an architectural flow diagram that 
25 illustrates general functional elements of the network code constructor. 

2.2 System Architecture: 

The processes summarized above are illustrated by the general system 
30 diagram of FIG. 2. In particular, the system diagram of FIG. 2 illustrates 
interrelationships between program modules for implementing a "multicast 

-13- 



network code constructor" for facilitating network based coding in a multicast 
environment by computing efficient multicast codes for optimizing network flows. 
As shown in FIG. 2, the multicast network is presented as a directed acyclic 
graph (DAG) 200 having zero-delay edges. However, as described in detail 
5 below, the ideas described herein are easily extensible to robust network codes, 
more general encoding scenarios, and networks with delays. 

As illustrated by FIG. 2, a network code constructor for computing efficient 
multicast codes for optimizing network flows begins by inputting known 

10 parameters of the network 200 into a network parameter input module 210. 

These parameters include information describing the sender, S, each receiver, R\ 
and R 2 , each internal node in the network, and each edge in the network along 
with capacities on each edge. The network parameter input module 210 then 
provides these parameters to a network reduction module 220. The network 

15 reduction then reduces the network to a form that is more conducive to analysis. 
In particular, the network Q 200 is reduced to a network whose edges have unit 
capacities by replacing each edge having a capacity c with c edges having unit 
capacity. 

20 Following the reduction of the network, a network capacity module 230 

uses a conventional max-flow algorithm to determine a maximum flow C, from the 
sender to each receiver R t . The multicast capacity C of the network is the 
minimum value of these maximum flows, i.e., C = min/ C f . Given a transmission 
rate R less than or equal to this multicast capacity C, the next step is to compute 

25 encoding vectors and decoding matrices for transmitting symbols at rate R from 
the sender the receivers across the network, and then decoding those symbols at 
each receiver. These codes are computed by a network coding module 240 
given the network parameters and the transmission rate R. Once computed, 
these codes 250 are stored for use in transmitting, coding and decoding symbols 

30 on the network 200. 
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3.0 Operation Overview: 



The above-described program modules are employed in a multicast 
network code constructor for facilitating network based coding in a multicast 

5 environment by computing efficient codes for optimizing network flows. The 
following sections provide a detailed operational discussion of exemplary 
methods for implementing the aforementioned program modules. In particular, 
the following description provides a discussion of initial definitions and 
background for explaining the network code constructor. Following these 

10 definitions and background, several theorems underlying the implementation of 
the network code constructor are provided for purposes of explanation. Finally, 
an overall summary of the system is presented in terms of a general system flow 
diagram. 

15 3.1 Definitions: 

The following discussion provides a description of definitions and 
assumptions used in implementing the network code constructor with respect to 
an acyclic network. As noted above, the descriptions of the network code 
20 constructor described herein are explained in the context of an analysis of 
multicast networks over directed acyclic graphs with zero-delay edges. A 
network is called cyclic if it contains directed cycles; otherwise it is called acyclic. 

In particular, in the following discussion, a network is represented by a 
25 directed acyclic graph (DAG) & = (V,E) with a vertex set V and edge set E . All 
edges are directed, and multiple edges are allowed between two vertices. The 
vertex which is at the head of an edge e shall be denoted by head{e) , and the 
vertex at the tail by tail(e) . The in-degree of a vertex v is defined as the number 
of edges e 1 such that head(e') = v , and is denoted by T(v) . The set of edges 
30 incoming to a vertex v is the set of all edges e such that head(e) = v , and the 
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incoming edges are denoted by {^(v),«f (v)e£ v) (v)}. Each edge eeE is 

assumed to have integer edge capacity and be delay-free. Note that as is well 
understood by those skilled in the art, by choosing a sufficiently large unit of time, 
any network can be approximated to an arbitrary degree of accuracy by a 
5 network with edges having integer capacities. 

For purposes of explanation, the following discussion describes the 
network code constructor in multicast networks, i.e., networks that have a single 
source vertex s and N receiver vertices or nodes 71 = {r v r 2i ...,r N } . Statistically, 

10 at every unit of time / , the source vertex s produces R binary random variables 
{bij} R M such that {6 /y .} are independently and indistinguishably distributed (i.i.d.) 
with Bernoulli (1/2) for all i e {0,1,...} , and j e {l,2,...,i?} . Further, for every group 
of m time units {im,im + l,... } im+m-\} the network codes described herein 
group together the m bits {4^,* MJ v-m&*^ u } to obtain R symbols {^,}* =1 in 

15 a finite field, IF 2 m , where m is some integer whose value will be fixed later. In 

addition, in order to simplify the notation in the following discussion, each edge is 
broken up into parallel edges with the unit capacity of 1 bit per unit time. 
However, it should be appreciated by those skilled in the art that the network 
coding methods described herein are equally applicable to any capacity of U bits 
20 per unit time. Therefore over m time units, each edge of the network has the 
capacity to transmit m bits of information, or equivalently, 1 symbol from the finite 
field IF 2 m . 

Therefore, in view of the preceding discussion, the network coding 
25 problem can be described as involving the transmission of the {X i} } R M symbols 
in m consecutive time units {0,1,..., m-l} over the edges of the network, to 
replicate {X 0J }*= X symbols at vertex v for every vertex veil. The identical 

process is then used for each successive block of m time units. Therefore, in 
order to again simplify the notation for purposes of discussion, the time index, i, 
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will be dropped from the notation such that symbols are denoted as {Xj} R M . 

Further, again for purposes of explanation, it is assumed that all edges are 
synchronized with respect to symbol timing, and that edges {<?_,, e_ 2 ,. are 
appended to the network such that head(e_j) = s for all j e {1,2,. . 

5 

Thus, a symbol Y(e_j) transmitted over edge e_ } is denoted by the 
symbol X. . Further, the symbol transmitted over any other edge e in the 
network, denoted by Y(e) , is a function only of the symbols on edges incoming to 
tail(e). Therefore, for each edge e , a Y(tail(e)) -length combination vector is 
10 defined as fi(e) = (A(e),/5 2 (e),...,^ r(to//(e)) (e)) to be associated with e . These 
combination vectors, (3(e) , represent encoding vectors for each edge of the 
network. The network g is said to be a JF 2 m -linear network if for all edges 
e e E , the symbol Y(e) on an edge e e E can be written as: 

15 Y(e) = Eqn. 1 

for some encoding vector (3(e) whose entries are elements of JF 2 m , where 
£?;"(v),...,4 n (V )(v) denote the edges incoming to a vertex v. These vectors (3(e) 

can be taken to be time-varying or time-invariant, and the network is accordingly 
20 either a "time-varying" or "time-invariant network. The following discussion 
concentrates on time invariant networks for purposes of explanation; however, 
the multicast network code constructor is clearly not limited to time-invariant 
networks, as it can easily be applied to time-varying networks. 

25 Furthermore, each edge e is associated with an R -length representation 

vector 77(e) , such that 
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77(e) = /3(e) 



77(er(tail( e ))) 
r7(e; n (tail(e))) 



Eqn. 2 



T(UO(e)) 



(tail(e)))j 



In addition, for each j e {1,2,..., , rj(e„j) is initialized as an R -length 
vector with 1 in the j th position and 0 everywhere else. Consequently, because 



5 of the directed acyclic nature of g , this initialization makes rj(e) well defined for 
each e e Af . Further, by using this definition of 77(e) , it can be seen that if the 
symbol Y(e_j) being transmitted over edge e_ } equals the symbol X. , then the 
symbol being transmitted over edge e equals 



Thus, the representation vector 77(e) on edge e shows how the symbol Y(e) on 
edge e is represented in terms of the original source symbols. 



Further, for each receiver r n R distinct T^) -length decoding vectors 

are defined to be associated with r r The entries of e(r i9 j) are 

elements of IF 2 m . Again, for purposes of explanation, the following discussion is 

restricted to the case where decoding operations also involve only linear 
20 combinations of the incoming symbols, as illustrated by Equation (4). However, 
as should be appreciated by those skilled in the art, the network code constructor 
described herein is also applicable to non-linear combinations of the incoming 
symbols. 



10 



Y(e) = V (e)(X v X 2 ^X R ) T 



Eqn. 3 



15 



25 



X(r»j) = e(r i9 j) (y^ft), W(Ov..^< tol7(e)) (^)) r 



Eqn. 4 
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For a given network & , the multicast network coding problem is defined 
as a 4 -tuple (<7,s,ft,R). Further, it is desired to choose appropriate P(e) and 

e(r iy j) so that the symbols can be decoded perfectly, i.e., X(r n j) = X } , for each 

receiver /;. € 71 and for each block of m consecutive time units. If such 0(e) and 

5 e(r„j) exist in a suitably chosen field JF 2 m , then the network & is said to be 

u JF 2 m -linearly solvable with multicast rate R", and the (5(e) and e(r n j) comprise 

the solution. 

Next, for the communication network 0 ,acut between two vertices v and 
10 V is defined as a partition of the vertex set V into two sets S and S c = V - S 
such that v e S and v'eS c . The value V(S) of the cut is defined to be equal to 
the number of edges going from S to S c . 

3.2 Theorems: 

15 

The following paragraphs describe several theorems that are utilized in 
formulating the derivation of algebraic network codes as provided by the network 
code constructor described herein. 

20 3.2.1 Theorem 1 - Min-Cut Max-Flow: 

The well known conventional "Min-Cut Max-Flow theorem" is described by 
L. R. Ford, Jr., and D. R. Fulkerson in "Maximal Flow Through a Network", 
Canadian Journal of Mathematics, 8, pp. 99-404, 1956. This Min-Cut Max-Flow 
25 theorem guarantees that for N = 1 , (that is, if there is only one receiver r { in the 

network) the network problem is solvable if and only if R is less than or equal to 
C ( . , the minimum value of all cuts between s and /;. . In particular, the Min-Cut 

Max-Flow theorem in network coding states that the network problem (<7,s,r„R) 

is solvable if and only if R < C ( . = min^^ V(S) . 
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A number of low time-complexity schemes have been devised based on 
the Min-Cut Max-Flow theorem for finding solutions to network coding problems. 
In general, such schemes are commonly referred to as "maximum flow 
algorithms." A number of these maximum flow schemes are described by: R. K. 
5 Ahuja, T. I. Magnanti, and J. B. Orlin, in "Some Recent Advances in Network 
Flows", SIAM Review, vol. 33, pp. 175-219, 1991. The network code constructor 
described herein makes use of such a maximum flow algorithm in a pre- 
processing stage. 

10 Fort example, in a conventional maximum flow algorithm, a directed path 

P = {e,,e 2 ,...,e„} between s and r t comprises an ordered set of edges from E 

such that tail(e x ) = s , head(e n ) = r t , and head^) = tail(e M ) . A set of C edge- 
disjoint directed paths P. between s and r t comprises C directed paths {P i; } 
between s and r. t such that no two paths share any edge in common, although 
15 sharing of vertices is allowed. Conventional maximum flow algorithms provide a 
means of finding C t edge-disjoint directed paths between s and r r In other 
words, these conventional maximum flow algorithms are used to find a solution 
for the network in which every coefficient in all the f3(e) and all the e^j) are 

either 0 or 1 . 

20 

3.2.2 Theorem 2 - Network Multicast Theorem: 

For the multicast network coding problem, let C be the minimum over all i 
25 and all cuts between s and r., that is, 

C = min min V(S) Eqn. 5 

In accordance with the conventional network multicast theorem described 
30 by R. Ahlswede, N. Cai, S.-Y. R. Li, and R. W. Yeung, in "Network Information 
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Flow", IEEE Transactions on Information Theory, IT-46, pp. 1204-1216, 2000, the 
multicast network coding problem (g,s,ft,R) is solvable if and only if R < C. 
One direction of the network multicast theorem (solvability of the network implies 
R < C ) is a direct consequence of the conventional Min-Cut Max-Flow theorem. 

5 The other direction requires more work. Proofs of this conventional theorem are 
provided by R. Ahlswede, N. Cai, S.-Y. R. Li, and R. W. Yeung, "Network 
Information Flow", IEEE Transactions on Information Theory, IT-46, pp. 1204- 
1216, 2000; and by R. Koetter, and M. Medard, in "An Algebraic Approach to 
Network Coding", presented at INFOCOM 2002 

10 (http://www.mit.edu/people/medard/pub.htm). 

3.2.3 Theorem 3 - Network Multicast Algorithm: 

15 

A conventional network multicast algorithm is described by R. Koetter, and 
M. Medard, in "An Algebraic Approach to Network Coding", presented at 
INFOCOM 2002 (http://www.mit.edu/people/medard/pub.htm). In that 
description, Koetter, and Medard explain that the multicast network coding 
20 problem {(/,s,7l,R) , with R < C, is "IF 2 m -linearly solvable" with multicast rate R 

if ro = [log 2 (iVC + l)l. 

While the network multicast algorithm described by Koetter and Medard is 
useful, a reduction in the field size, as accomplished by the network code 
25 constructor, as described herein serves to increase network throughput without 
flooding the network. 

3.2.4 Theorem 4 - Improved Network Multicast Algorithm: 

30 

The network code constructor described herein provides an improved 
network multicast algorithm that provides for multicast network coding in which 
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arithmetic operations can occur in any finite field with more than AT-1 elements, 
where N represents the number of receivers in the network. 

In particular, as described in detail below, the network code constructor 
5 described herein demonstrates that the network coding problem (£,s,7l,i?) , with 

R<C (Theorems 1 and 2), is "IF 2 m -linearly solvable" with multicast rate R if 

m = flog 2 N] . This is a significantly smaller field size than has been 

demonstrated by any other conventional method. Clearly the advantages of a 
smaller field size include decreased computational requirements, fewer symbols, 
10 and greater network throughput. 

Thus, the remainder of this description of the network code constructor 
assumes R<C (Theorems 1 and 2) and that m equals [log 2 N}. In particular, 

15 Let the flow to r. , 7> = be a set consisting of R edge-disjoint directed 

paths P. j between s and . A network flow Af(R) between s and TZ is defined 

as the union of the edges in the flows to each of the receivers, i.e., 

Af(R) = {eeE: 3P ijy e e P itJ } . The directed acyclic nature of Q induces a total 

order < on the edges in Af(R) such that for any edges e and e 1 , 

20 head(e) = tail{e f ) implies e < e . 

3.3 Multicast algorithm for Zero Delay Directed Acyclic Graphs: 

In view of the definitions and theorems provided above, the processes 
25 utilized by the network code constructor can now be described. In particular, the 
following paragraphs describe the multicast network code design computed by 
the network code constructor in terms of algorithm inputs, algorithm outputs, 
preprocessing, calculation of encoding vectors (3(e) for every edge e of the 

multicast network, and calculation of decoding vectors e(r n j) for decoding 
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symbols multicast from the sender to each receiver. In general, the encoding 
vectors comprise a multicast network encoder, and the decoding vectors 
comprise a multicast network decoder. The multicast network encoder and the 
multicast network decoder are collectively referred to as a "multicast network 
5 coder." Further, the term "multicast network code" is used to refer to the 
collection of encoding vectors on each edge, and decoding matrices at each 
receiver. 

3.3.1 Code Design: 

10 

The input needed by the network code constructor is simply the 
aforementioned 4 -tuple (Q,s,Tl,R), namely a DAG G, labeled with edge 
capacities, the sender node s , receiver nodes TZ , and the multicast rate R. As 
described in detail below, given this 4 -tuple input, the network code constructor 
15 outputs encoding vectors P(e) for every edge e in the network, along with 

decoding vectors {e(r X9 j)}f ml for use in multicasting symbols from the sender and 
decoding those symbols at each receiver. 

As illustrated by FIG. 3, the network ff\s a DAG having known nodes and 
20 edges 300, with integer capacities on each edge. The first step in determining 
codes for the network is to reduce the network to a form that is more conducive 
to analysis. In particular, the network (/ is reduced to a network whose edges 
have unit capacities by replacing each edge having a capacity c with c edges 
having unit capacity 310, and output 320 for use in determining network codes. 

25 

Next, any conventional low time-complexity maximum flow algorithm, as 
noted above in Section 3.2.1 , is run N times (once for each receiver node) to 
obtain a network flow M(R) between the sender s and the set of receivers/? . In 
particular, as illustrated by FIG. 4, the network f 400 having c edges with unit 
30 capacity output as illustrated by FIG. 3, is used as in input for determining the 
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multicast capacity of the network. Specifically, for each receiver node r { , a 
conventional max-flow algorithm is used to find the max-flow C ( . from the sender 
node s , to r t 420. After the max-flow C,. each receiver node r t is computed 420, 
i is incremented for selecting the next receiver 430. 

5 

This process (steps 410 through 430) repeats until all Af receiver nodes 
have been processed. At that time, the multicast capacity C of the network 
between s and TZ is determined by simply identifying the minimum of the 
C, values computed for each of the N receiver nodes 440. This minimum of the 

10 C ( . values is then output 450 as the multicast capacity C of the network between 

s andft. 

3.3.2 Computing Encoding Vectors: 

15 Given the network 0 with unit capacity edges, and the multicast capacity 

C of the network between s and 72. , the next step is to choose a transmission 
rate R£C, and to compute a flow 7? consisting of R edge-disjoint paths 
Pi j e V, from s to each r. t in K . Note that the resulting network flow N(R) will 
restrict the network code to utilize only edges in N{R) , rather than possibly all 

20 edges in E . The next step is then to determine encoding vectors, p(e t ) , and 
decoding vectors, ri(e t ) , for the multicast network. In general, this is 
accomplished by finding (3(e t ) and r){e t ) for te{\,...,\J\f(R)\}. At each step /, a 
frontier edge set for receiver r { at step / is defined as an ordered subset of 
R edges from 7? such that the j ^ edge is the edge e t , in the ;' th path P ( j e 7? with 

25 the largest t' less than or equal to /. A frontier edge set matrix for receiver r. at 
step t is an RxR matrix whose rows are the representation vectors r}(e,) 
for the edges in ^ (<) . 
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Therefore for each step t , there are k (not necessarily distinct) edges 
{e lt ,e l2 ,...,e lk } in the frontier edge sets Jj'"'.*,.., , such that at step / 
the network code constructor replaces the edges {e f ,e, 2 ,...,e, t } in the frontier edge 
sets in which they appeared with e, , to obtain the updated frontier edge sets. 
5 The network code constructor then calculates for e, an encoding vector (3(e t ) , 
and Equation (2) is used to obtain the vector ri(e,) . The frontier edge set matrices 
{^'"'^j'^...,^' -0 } are also updated by replacing the vectors 
{v(e li ),ri(e lj ),...,r)(e tk )} with j](e,). In particular, the encoding vector /?(e,) is 
chosen to be any combination vector so that 

10 

Vi e {1,2,..., W}, rank(T^) = R Eqn. 6 

The step-counter / is then incremented by 1 and this procedure repeats until 
t =| M(R) | . Once this procedure has been completed, each frontier edge set 

15 consists only of edges e such that head(e) = r t , and consequently, all 

encoding vectors, /3(e) , have been determined. 

Note that for purposes of explanation, a proof that P(e,)'s can be chosen 

to satisfy Equation (6) is provided in the following paragraphs, along with a more 
20 detailed explanation of the computation of the encoding vectors, 0(e) . 

In particular, the existence of/?(e)'s that satisfy Equation (6) is 
demonstrated in accordance with the following discussion. Specifically, for all 
i e {l,2,...,k} lei Mi be arbitrary non-singular nxn matrices over the finite field 

25 W q . Further, for all i e {1,2,...,£} , let v . be some row vector of M. . Let £ be 
any subspace of (TF Q ) n containing the span of {v,,...,v t } . Finally, for all 
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i e {1,2,..., k) , let S/ be the linear span of all the row vectors of M, except v , , i.e. 

Sj = span{rows(M .) — v ( .} . 

Given the above definitions, it is shown below that if k < q , then there 
5 exists v e £ such that for all i e {1,2,...,/:} , S t © v = (IF 9 )" . Specifically, it should 
be noted that C-\J M S, equals £-u* =1 (£nS ( ) , and that by using a counting 
argument on the set £-u* =1 (£n5' ( .) the required result is obtained. 

In particular, let d(£) be the dimension of the vector space £ . Therefore, 
10 the number of vectors in £ equals q d(C) . For each i e {1,2,..., A:} , the dimension of 
the vector space (£n5,) is strictly less than d(C) , since £ contains at least one 
vector, v ( . , which is not contained in 5 ; . Therefore, for each j e {1,2,..., A:} the 
number of vectors in £nS ( . equals at most q d{C) ~ l . In addition, each £nS ( . 
contains the zero vector, and therefore the total number of vectors in U* =1 (£nS f ) 
15 is no more than k{q d(C) ~ x - 1) + 1 . Since k < q , this quantity is strictly less than 
the number of vectors in £ , q i(C) . Therefore | £-u* =1 5, |> 0 , which implies the 
existence of v e £ such that v £ S, for any i e {1,2,...,&} . This observation 
provides the desired result. In fact, the converse is also true, in the sense that it 
is possible to construct # + 1 subspaces S,. such that U^S f = (IF^)" , and 

20 therefore there is no vector v with the desired properties. 

Further, by extension, a corollary to the above proof is provided to show 
that for all eeM(R) , there exists a p(e) such that Equation (6) is satisfied. In 

particular, set n = R and q = 2 m . Let edges {e, ,e (2 ,...,e, t } be the edges in the 

25 frontier edge sets {^ i ( '" ,) ,^ ( '" 1) ,...,^ t ( '" 1) } , which are replaced at step t by e,. 

Then M j = Fy~ i) for j e {i v i 2 ,...,i k } . Finally, let £ be the span of 

{iK«f(totfte)))>E? } . 
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Given the above proofs, and general description of the computation of the 
encoding vectors, /3(e t ) , an algorithm on the order 0((R + rfr) is provided for 
finding /?(e,)'s for any edge e t ^M(R). Specifically, let v., S,, £, d(C), and S, 
5 be as defined in the preceding paragraphs. Then, given inputs of row vectors 
{v,}/ = i and (R-l)*R matrices {5,}f =1 , outputs of row vectors v are computed by 
the multicast network encoder such that v e £ , v £ S, for any i € {1,2,..., A:} , along 
with row vector 0 = (0 v /3 2 ,...,p k ) such that v = J^* =] #.v ( . . 

10 In particular, let L be a d(C)xR matrix whose rows form a basis for £ . 

For each i e {1,2,..., A:} , it is possible (by row operations on L and the matrix 5, ) to 
obtain the (d(C)-\)xd(C) matrix B, whose rows form a basis for £nS ( in the 
coordinate system given by the rows of L . For each /' , a d(C) -length column 
vector Zi is identified, called the zero vector of £n 5 f in (W q ) n , such that for any 

15 vector y' in the span of the rows of B i , y'. z f = 0 . To obtain such a vector, row 
operations on are performed on B, until it contains the (</(£) -l)x(rf(£)-l) 
identity sub-matrix and a (</(£) -1) -length column vector z' . z/ then equals the 
row vector obtained by adjoining the element 1 to -z'. , i.e., z , = (-z' r >l) • 

20 A d(£) -length row vector y is then identified such that 

y. z, * 0 for any i € {1,2,..., A:}, Eqn. 7 

and an /? -length vector v = yL is identified in £ such that v £ £nS f for any 
25 / € 1,2,..., A: . Given such a v , 0 t is easily computed by a process of Gaussian 
elimination. A vector y = (y v y 2 ,-,y d (C)) satisfying Equation (7) can be obtained 
via a "greedy algorithm" as described below for use in computing /?, . 
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3.3.3 Greedy Algorithm: 

A "greedy algorithm" for computing the vector y = (^j 2r ,^ (£) ) satisfying 
Equation (7) operates by first denoting a kxd(C) matrix whose i th row vector is 

5 zf by Z . Without loss of generality, let the first rank(Z) rows of Z be linearly 
independent, and denote this rank(Z)xd(C) sub-matrix by Z ; . Further, denote 
the matrix consisting of the remaining rows of Z by Z D . Therefore Z D can be 
written as the matrix product TZ l , for some {k-rank(Z))xrank{Z) matrix T . 
The greedy algorithm described herein computes a rank(Z) -length column 

10 vector c such that any vector y satisfying Z 7 y r = c satisfies Equation (7). 

In the first step of the greedy algorithm some arbitrary value c x ^ 0 (such 
as, for example, c x = 1 ) is chosen. Now there are two possibilities - either the 
rank(Z) = 1 , or it is greater than 1 (it cannot be zero, since all the row vectors of 
15 Z are non-zero). If rank(Z) = 1 , c is assigned the value (1) . If rank(Z) > 1 , 
more work is needed to calculate c . For notational convenience the / th row 
vector of the matrix T is denoted by (T) i , and the (ij)** element by (T) u . 

Consider all row vectors (T). of T which have non-zero elements only in the first 
ctr positions, and denote them by the superscript ctr , as (T)f . where ctr is a 
20 counter that is initialized to 2. The greedy algorithm proceeds inductively in the 
variable ctr and at each step of ctr computes a constant, c ctr e JF 2m , such that 

(Z 7 ) c/r y r =c c/r ^0 Eqn.8 

25 is true. To choose this value c c , it is noted that 
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ctr-\ 

# <Z#)*y * -(n- c V(E(r), y .(Z / ).)y r Eqn.9 

cfr-1 

where the rf,. c/r are some constants in TF 2 m - But since rank(Z) > 1 , there are at 
most £ — 2 row vectors in Z^ . However, since k < q = 2 m , there are at most 
5 2 m - 2 different values of d i ctr for a fixed value of ctr . Therefore there exists at 

least one c ctr such that (Z f ) ctr y T = c cfr does not contradict Equation (8). This 
process is continued until ctr = rank(Z) , at which point under-determined linear 
equations Z 7 y r = (c v c 2 ,...,c rank(Z) ) T = c. are produced. Therefore choosing any 

pseudo-inverse Z7 1 of Z ; and evaluating Z/'c 7, gives a y r such that the column 
10 vector Zy T has no zero elements. 

3.3.4 Computing Decoding Vectors: 

After computing the encoding vectors, /3(e) , as described above, the 
15 decoding vectors, e(r l9 j) , are then computed as described below. In particular, 
for all ie {1,2,..., N}, ye{l,2,...,i?}, e(r 99 j) is chosen as the j th row of the matrix 
(pwmyi This definition is well-defined since by assumption, as discussed 
above in Section 3.3.1 , the frontier edge set matrices are invertible. Therefore, 
for any receiver r { : 
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where the equalities follow from the definitions of the decoding and 
representation vectors described herein. 

3.4 System Operation: 

5 

The processes described above are illustrated in the sequence of system 
flow diagrams provided as FIG. 5 through FIG. 10. This processes depicted in 
these figures are provided for purposes of illustration, and are intended merely to 
illustrate one way of implementing the network code constructor. Clearly, in light 
10 of the description provided herein, there are a number of ways for computing 
codes for a multicast network, and the network code constructor is not intended 
to be limited to the processes illustrated by FIG. 5 through FIG. 10. 

3.4.1 Computing the Network Code: 

15 

As illustrated by FIG. 5, in computing codes for multicast networks, the 
first step is to provide several inputs that define the multicast network to the 
network code constructor. As illustrated by FIG. 5, these inputs 500 to the 
network code constructor include a directed acyclic graph 0 , which has been 

20 reduced to unit capacity edges (see FIG. 3 and associated discussion). In 
addition, the inputs to the network multicast coder also include a transmission 
rate R that is less than or equal to the maximum multicast capacity C of the 
multicast network (see FIG. 4 and associated discussion). Finally, the inputs to 
the network multicast coder include a known source node s, and known receiver 

25 nodes n of the multicast network. 

Given these inputs, the first step is to initialize the network code 
constructor as illustrated by the steps shown in FIG. 5. In particular, once the 
aforementioned inputs have been provided, the first step is to create a new node 
30 denoted by s f , and new edges denoted by e_ l9 ... 9 e_ R from s ' to the known 
source node s 505. Next, for each edge 510, e_ } , an TMength representation 
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vector 77(e_ y ) (where, as described above, R represents the transmission rate) is 
set equal to the i?-dimensional row vector, with 1 in the / h , and 0 otherwise 515. 
After rj(e_j) is set for each edge, e_j J is incremented for selecting the next edge 

520. In other words, for each je {1,2,.. .,/?}, r)(e_j) is initialized as an R -length 
5 vector with 1 in the y th position and 0 everywhere else. As noted above, because 
of the directed acyclic nature of ^ , this initialization makes rj(e) well defined for 

each eeAf . This process (steps 510 through 520) repeats until all j edges for s 1 
have been processed. 

10 Next, for each receiver 525, r { , a conventional flow algorithm is used to 

compute a flow of magnitude R (i.e., R edge-disjoint paths) from s 1 to /;. 530. 

After the flow is computed 530 for each receiver r { , i is incremented for selecting 

the next receiver 535. This process (steps 525 through 535) repeats until all / 
receivers have been processed. Next, any edges not in the flow are simply 
15 discarded 540. In other words, any node having excess capacity will have edges 
that are not utilized. These unutilized edges are discarded 540 for purposes of 
computing multicast network codes. Finally, the remaining edges, e p ...,e r , are 

ordered topological^ from s 545. This step completes the initialization stage of 
the network code constructor. 

20 

As illustrated by FIG. 6, the process continues (see box A (550) in FIG. 5 
and FIG. 6) by computing encoding vectors (3(e t ) for the network code 
constructor in an "encoder design stage." First, the total number of edges T 
representing the edges, e p ...,e r , computed in the initialization stage, is used as a 

25 counter 600 for computing the encoding vectors p(e t ) for each edge in a nested 
loop. In particular, for each t = {1, 2, 7} 600, a counter k is initialized to zero 
602. This counter k serves to keep track of the number of receivers whose flows 
use edge e t . Next, for each receiver r { 605, the processes for computing the 
aforementioned frontier edge set for each receiver are initialized 610 by defining 
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Pjj as the j th path (/ =1, R) in the flow to n\ defining e u as the edge e v in P u with 
the largest t' = t\ and defining = {e iv e iR ) as the "frontier edge set" for 
receiver r, at step /. 

5 Next, a determination 615 is made as to whether the current e t is in the 

frontier edge set, T x . (,) , i.e., e t = e ip for some j*. If is in then, as 
illustrated in box 620, the counter k is incremented by 1 to indicate that the 
current edge e t is used by the current receiver r h Further, if e t is in f then 

y k = n(e f ) , where e 1 is the predecessor of eij* in path P ir \ further, 

io S k = M^X.^ry^.J^^^^^^ , where {vJuSJs a linearly 

independent set of vectors. Then, whether or not e t is in , i is incremented 

for selecting the next receiver 625. This process (steps 605 through 625) 
repeats until all i receivers for the current edge / have been processed. 

15 Next, as illustrated in box 630, a linear combination of v p ..., is chosen 

630 such that rj is not in the span of any S X9 ...,S k , where {77} US,. is a linearly 
independent set of vectors for/ = . Note that this process for choosing a 
linear combination of v p ..., 630 is illustrated in greater detail with respect to 
FIG. 8, as described below. Next, as illustrated in box 635, r?(e,) is set equal 

20 to^, and P(e t ) is set equal to(/3 p ...,/^) . Then, after setting rj(e t ) and (3(e t ) 635, 

t is incremented for selecting the next edge 640. This process (steps 600 through 
640) repeats until all T edges have been processed. This step completes the 
encoder design stage of the network code constructor. 

25 As illustrated by FIG. 7, the process continues (see box B (645) in FIG. 6 

and FIG. 7) in a "decoder design stage" by computing a matrix of decoding 
vectors e^) at each receiver r x and outputting these decoding vectors along with 

a reduced directed acyclic graph (representing the multicast network) having an 
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encoding vector /?(<?) on each edge. First, for each receiver^. 700, e(r)) is 
determined by computing 710 a corresponding inverse matrix of decoding 
vectors sty) as illustrated by Equation 1 1 which follows from Equation 10 as 

described above in Section 3.3.2: 



Eqn. 11 



This process loops while incrementing 720 the current receiver r { until all 

receivers have been processed. Once all of the receivers have been processed 
10 (700 through 720), then the network code constructor outputs 730 the decoding 
vectors e(r)) along with a reduced directed acyclic graph (representing the 

multicast network) having an encoding vector (3(e) on each edge. 
3.4.2 Choosing a Linear Combination of Vectors: 

15 

As noted above with respect to Box 630 of FIG. 6, FIG. 8 illustrates the 
process for choosing a linear combination of v,,..., \ k such that 77 is not in the 

span of any S v ...,S k , where {77} u -S,. is a linearly independent set of vectors for 

1 = l,...,ft . 

20 

In particular, as illustrated by FIG. 8, the linear combination of v p ...,v 4 is 
chosen by first inputting the tf-dimensional row vector and the set of R-l R- 
dimensional vectors S t for i = l,... f * 800. Next, denoting by L the vector space 
spanned by v. for i = , the set of vectors v f is computed for i = , 
25 spanning Z,, where k 1 is less than or equal to k 81 0. In the next step, for 
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j = 1,...,* a vector Z ; is computed by Gaussian elimination in L such that for any 
vector y in Sj , y.Sj = 0 (Box 820). 

Once the vector Z, has been computed 820, the next step is to use this 
5 information to find a vector v in L such that \.Z. * 0 for j = l,...,k 830. Further, 

in finding this vector v in L it is assumed that vectors Z y . , j = 1 k" , are linearly 

independent and that all other Z. , j = k" + 1,...,* , can be written as linear 
combinations of these vectors. In the next step (Box 840), v is determined by 
finding a linear combination of vectors for j = l,...,k" by solving a system of 

10 linear equations denoted by v.Z 7 . = c,. for j = \,...,k" . Note that the solution to 

this equation may not be unique, and that one of the solutions may be calculated 
by Gaussian elimination. Finally, once this system of linear equations has been 
solved 840, an output of k coefficients /?,,. and the i?-dimensional vector 

rj = ^2 tJ t 0,y, is provided. This output is then provided to Box 625 of FIG. 6 as 
15 noted above. 

In addition, the step illustrated by Box 840 (finding a linear combination of 
vectors c y for j = \,...,k" ) is further described with respect to FIG. 9. In 

particular, as illustrated by FIG. 9, finding a linear combination of vectors c y for 

20 j = \,...,k" begins by inputting the linearly independent vectors Z,., i = \,...,k" , 

and Zj, j = k" + 1,..., k dependent vectors on the first vectors 900. 

Next, for each vector, j = k" + \,...,k , a dependence number is defined as 
the largest i such that Z y .is dependent upon Z ( . 910. Then, as illustrated by Box 

25 920, for each i = \,...,k" , c. is fixed, thereby fixing all of the v.Z y for each j having 
a dependence number i. Further, a non-zero c ( . is then chosen so that none of 
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the \.Zj fixes to zero for each j having a dependence number L This choice is 
easily made by simply considering k—k" different values of c, to identify the 
proper selection of each c. . Having made this selection 920, c,,...,c A is output 
930. As noted above, this output is provided to Box 840 of FIG. 8. 

5 

3.5 Multicasting Symbols from the Sender to Every Receiver: 

Once the aforementioned encoding vector and decoder matrix have been 
computed for the multicast network, then they can be used for transmitting 
10 symbols across the network from the sender to be decoded by each receiver. 
This process is illustrated by FIG. 10. 

" In general, as illustrated by FIG. 10, the multicast network coder (i.e., the 
encoder and decoder) allows for the efficient multicasting of a set of symbols 
15 from the source to every receiver in a multicast network. The inputs used to 
produce this multicasting include a rate R < C , a set of symbols X v X R to be 
multicast, a directed acyclic graph = (V,E) of nodes and (unit-capacity) edges, 
a source node s in V, N receiver nodes r v r N in V, an encoding vector (3(e) on 
each edge e, and a decoding matrix e(^) at each receiver n (Box 1000). Given 

20 these inputs 1000, the next step is to create a new node denoted by s f , and new 
edges denoted by e_ l9 ...,e_ R from s ' to the known source node s 1005. Next, for 
each new edge 1010, an input symbol Y(e_j) is set equal to X r This 

process is repeated for each new edge 1020 until all new edges have been 
processed. 

25 

Then, in topological order from s for each remaining edge e 1025, Y(e) is 
computed on edge e 1030 by defining v = tail(e) as denoting the node that e is 
leaving, and defining e X9 ...,e k as denoting the edges entering v. Given these 
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definitions, the computation of Y(e) 1030 is accomplished by evaluating the 
following expression, derived from Equation (1): 



r<«) = fe(«) - 



Eqn. 12 



This process (1025 through 1035) is repeated for each remaining edge 1035 until 
all of the remaining edges have been processed, thereby encoding the input 
symbol Y(e_j) being multicast by the sender. 



10 Next, for each receiver, r, 1040, the symbols X v X R are computed 

1045 by evaluating the following expression, derived from Equation 4: 



Eqn. 13 



15 This process (1040 through 1045) is repeated for each receiver, r, 1050, until all 
of the receivers have been processed. 

Finally, having computed the symbols A",, X R at each receiver 1045 

those symbols are output 1055, thereby completing the multicast from the sender 
20 to all of the receivers by having decodes the symbol multicast from the sender to 
all of the receivers. 



In conclusion, it should again be noted that although the network code 
constructor described herein is presented in terms of a low complexity process 
25 for computing multicast codes for networks in an information theoretically optimal 
manner, the ideas described herein are easily extensible for obtaining algebraic 
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codes over small finite fields for more general networks, networks with delays 
and robust networks. 

The foregoing description of the invention has been presented for the 
5 purposes of illustration and description. It is not intended to be exhaustive or to 
limit the invention to the precise form disclosed. Many modifications and 
variations are possible in light of the above teaching. Further, it should be noted 
that any or all of the aforementioned alternate embodiments may be used in any 
combination desired to form additional hybrid embodiments of the audio 
10 challenger described herein. It is intended that the scope of the invention be 
limited not by this detailed description, but rather by the claims appended hereto. 
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